Realization of Levy flights as continuous processes 
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On the basis of multivariate Langevin processes we present a realization of Levy flights as a 
continuous process. For the simple case of a particle moving under the influence of friction and a 
velocity dependent stochastic force we explicitly derive the generalized Langevin equation and the 
corresponding generalized Fokker-Planck equation describing Levy flights. Our procedure is similar 
to the treatment of the Kramers-Fokker Planck equation in the Smoluchowski limit. The proposed 
approach forms a feasible way of tackling Levy flights in inhomogeneous media or systems with 
boundaries what is up to now a challenging problem. 
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It has become evident that Brownian random walks 
may be a too simple description of diffusion processes 
in complex systems like the motion of tracer particles in 
turbulent flows the diffusion of particles in random 
media 0, the motion of wandering albatrosses 01, hu- 
man travel behavior and spreading of epidemics J| or 
economic time series in finance [5|. A variable x cor- 
responding to such a process can frequently exhibit the 
dynamics described by the notion of superfast diffusion, 
where the characteristic value x of the variable x demon- 
strates scahng behavior oc t^/" with a <2 [1, 01- 

Brownian motion is described on the basis of Langevin 
equations or, in a statistical sense, by the Fokker-Planck 
equation (cf. [1, Q). A straightforward way to deal with 
anomalous diffusion is based on a generalization of the 
Langevin equations by replacing Gaussian white noise 
with Levy noise Recently, there has been a great 

deal of research about superfast diffusion. It includes, 
in particular, a rather general analysis of the Langevin 
equation with Levy noise (see, e.g., Ref. [nj) and the 
form of the corresponding Fokker-Planck equations (l^ . 
[l3j . description of anomalous diffusion with power law 
distributions of spatial and teraporal steps |,10l . [13] , Levy 
flights in heterogeneous media 15, 3, 13 ^^'^ external 
flelds [3, [l^ I fi'^^t passage time analysis and escaping 
problem for Levy flights 0, [U [H, [H [H [111 , as well 
as processing experimental data for detecting the Levy 
type behavior I26!|. Besides, it should be noted that the 
attempt to consider Levy flights in bounded systems (see, 
e.g., Ref. [27, 28 1 and review [2^ as well) has introduced 
the notion of Levy walks being a non-Markovian process 
because of the necessity to bound the walker velocity. 

The key point in constructing the mutually related 
pair of the stochastic Langevin equation and the non- 
local Fokker-Planck equation for superdiffusion is the 
Levy-Gnedenko central limit theorem [13, [3, [13] ■ For 
the superdiffusion processes it specifies the possible step 
distributions P{Ax) which are universal and actually in- 



dependent of the details in the system behavior at the 
microscopic level. In particular, for a symmetrical ho- 
mogeneous ID system superdiffusion can be regarded as 
a chain of steps {Ax} of duration St whose distribution 
function P{Ax) exhibits the following asymptotic behav- 
ior for \Ax\ 3> x{St) 



P{Ax) 



(1) 



In spite of the considerable success achieved in this field 
the theory of superdiffusion is far from being completed. 
For a given elementary step of any small duration it 
is impossible to single out some bounded domain that 
contains its initial Xi and terminal xt points with the 
probability practically equal to unity because the second 
moment ((xt — x^)^) diverges. This renders the descrip- 
tion of Levy flights in heterogeneous media or media with 
boundaries a challenging problem. Within the classical 
formulation the Levy flight is not a spatially continu- 
ous processes. As a consequence it is not possible to at- 
tribute local characteristics to Levy flights which might 
help to identify, e.g., the encounter time with boundaries. 
Bounding the particle velocity breaks the Levy as well as 
the Markov properties. 

The purpose of the present Letter is to formulate an 
approach to describing Levy flights and Levy walks using 
the notion of continuous Markovian trajectories. The key 
idea is to introduce the velocity as a second variable but 
staying with simple Gaussian noise. For a fixed time 
scale 6t we can recover the standard behavior of Levy 
flights. However, we have full locality in the sense that a 
trajectory can be determined with any desired resolution. 

In the general form the model proposed for considera- 
tion is reduced to the class of coupled governing equations 
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FIG. 1: Characteristic form of random walks described by 2D 
analogy of model (|3]). The used system parameters meet the 
Levy exponent a — 1.6. 



for the state vectors x = {xi} and v — {vi} 



dx 

m 

dw 



F(x,v), 

H(x,v) + G(x,v).|(t). 



(2a) 
(2b) 



Here the Langevin equation \2h\ is written in the Ito 
form, ^(t) = {S,i{t)} are the collection of mutually inde- 
pendent Gaussian white noise components, the "forces" 
F(x,v), and H(x, v) are given functions, and the matrix 
G(x, v) depending on the state variables specifies the in- 
tensity of Langevin "forces". In some sense we reduce 
superdiffusion to a normal diffusion process expanding 
the phase space where a new variable, particle velocity, 
is governed by the Langevin equation with normal but 
multiplicative noise. 

In this Letter we exemplify our procedure applying to 
the following ID system with two variables, the coordi- 
nate of random walker x and its current velocity w. 



dx 

-di=''' 

dv _ {a + 1) 

~ 2r 



(3a) 
(3b) 



Here r is a certain time scale, the intensity of the Lan- 
gevin random force is given by the function 



(4) 



with the parameter Va measuring the intensity of the ad- 
ditive component of Langevin force, ^(t) is white noise 



such that {£,{t)^{t')) = S{t — t'), and the parameter 
a € (1,2). The Langevin equation (j3bp is written in 
the Hanggi-Klimontovich form, which is indicated by the 
symbol *. The dynamics, resulting from a 2D version of 
these equations, is visualized in Fig. [TJ 

The corresponding forward Fokker-Planck equation for 
the distribution function V{x — xo,v,vo,t) reads 



dp 

dt 



2t dv 



dV 

g''{v)— + {a + l)vV 
ov 



(5) 



where the values xq and specify the initial position of 
the walker. 

The distribution of the walker velocities is deter- 
mined by the partial distribution function 



Pv{v,vo,t) := {d{v ~ v^)) 



(6) 



and, by virtue of ([5]), the stationary velocity distribution 
P^'(u) meets the equality 

whence we immediately get the expression 



(7) 



where r(. . .) is the Gamma-function. In addition, us- 
ing the Fokker-Planck equation for function ([6]) following 
directly from ([5]) we find the expressions 



{v{t)) = vo exp 



{a - 1) t 
2 T 

{2 -a)- 

T 



(8a) 

for Wo > Va (8b) 



characterizing actually the relaxation of the initial veloc- 
ity distribution to its steady state form. 

The found exponential decay of the first velocity mo- 
ment demonstrates the fact that the Levy walker "re- 
members" its velocity practically on time scales not ex- 
ceeding the value r. The exponential divergence of the 
second moment (j8bl) indicates that the system relaxes to 
the stationary distribution ([7]) on time scales t^r. So, 
in some sense, the spatial steps of duration about r are 
mutually independent. In other words, the value t sep- 
arates the time scales into two groups. On scales less 
than T the particle motion is strongly correlated and has 
to be considered using both the phase variables x and v. 
Thus, on a time scale 6t ^ t the particle displacements 
are mutually independent and the succeeding steps of the 
Levy walker form a Markovian chain, with the particle 
velocity playing the role of Levy noise. This scenario is 
exemplified in Fig. [2] for some realization of v{t) follow- 
ing from equation (j3bp . Levy flight events, i.e. the long- 
distance jumps of a Levy walker, are due to large spikes 
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FIG. 2: Characteristic form of the time pattern v{t) exhibited 
by the stochastic system (|3]). The individual windows depict 
the patterns on various scales. In simulation the parameter 
a — 1.6 was used. 



of the time pattern v(t) whose duration is about several 
T. More precisely, the long-distance displacement Ax of a 
walker during a certain time interval St is mainly caused 
by the velocity spike of maximal amplitude 1} attained 
during the given interval, i.e. Ax ^ "dr. For 6t ^ t the 
quantity is statistically uncorrelated during succeed- 
ing time intervals. 

Now we proceed in a two steps. First, we use this sim- 
ple physical picture to show via a combination of analyti- 
cal and numerical evidence that the distribution function 
Px{Ax,vo,t) := {5{x — xo — Ax)) indeed is of form ^ 
for t ^ T. Second, we strictly show that the correspond- 
ing generating function fulfills 



exp 



r(a)r f r 



(9) 



The latter expression is the standard generating function 
of Levy flights with exponent a and matches the distri- 
bution H]). 

If the spikes in Fig. [5] had the same shape and St t 
the normalized walker displacement Ax/-d would be a 
constant of the order of r (in the limit of large d where 
Aa; is largely determined by a single peak). Then Px{Ax) 
would directly follow from the distribution of maximum 
velocities. To proceed we, first, make use of the relation 
between the extremum statistics of Markovian processes 
and the first passage time distribution 3l|. Namely, the 
probability function <I>(i?,w,t) of the random variable d 
and the probability i^(i3,w,t) of passing the boundary 
V = ±z9 for the first time at moment t are related as 



H^,v,t) 



d_ 



dt'F{d,v,t') 



(10) 



Here v is the initial position of the Levy walker. Then 
analyzing the Laplace transform of the first passage time 
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FIG. 3: The ratio Ax I'd for individual steps vs the values of 
the random variable i9. In simulation a = 1.6 was used. 



distribution Fi^ [d, v, s) we will get the conclusion that the 
distribution function <i>('i9,w,t) of the velocity extrema 
is of the form (see the supplementary materials) 



_d_ 

W) 



(11) 



for V < 'd{t) and -d > ■i?(i). Here the quantity t?(t) = 
Va{t/T)^^°' is the velocity scale characterizing variations 
of the random value "d and the function (t>{Q possesses 
the asymptotics 
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(12) 



Via numerical simulation we have determined the dis- 
tribution of Ax I'd for given velocity extremum The 
first and second moment of this distribution is shown 
in Fig. [3l As expected the average value of Ax/d in- 
deed approaches a constant Ct (for a = 1.6 the value 
Cr ~ 1.6r). However, the finite variance shows that 
the velocity spikes have some distribution in their shape. 

and P{Ax) are 



Thus a priori the distributions 

not identical when replacing Ax by Cr'd- However, since 
the distribution of Axj-d for fixed 'd does not depend on 
d (for large d) one can directly write 



Px{Ax) cx j ded^q{e)d-^^+"'>S{e^ + Crd - Ax) 

cx / rfeg(e)[(e + c,)/Ax](i+"' oc Ax-(i+") , (13) 



where q[e) is the distribution of the random variable e := 
Axj-d — Cr- Thus, despite the variance in peak shapes 
the algebraic distribution of -d directly translates into an 
identical distribution for Aa;. 

We have performed a stick derivation of formula ^ in 
the following way (for details see the supplementary ma- 
terials available online). The appropriate Fokker-Planck 
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equation should be written for the full generating func- 
tion G{>c, fc, t) for system ([3]) which depends on two wave 
numbers, x as before and k related to the velocity vari- 
ations. Then the corresponding eigenvalue problem can 
be analyzed assuming the wave number >f to be a small 
parameter. It turns out that the perturbation caused by 
the K-teim is singular which affects essentially the mini- 
mal eigenvalue, making it dependent on x as Amin oc >ir". 
In this way expression Q is obtained. Furthermore the 
specific value of Cr equal to 



equation. For equilibrium systems the fluctuation dissi- 
pation theorem relates linear damping and purely addi- 
tive noise. The emergence of Levy flights, however, is 
related with the presence of multiplicative noise, and, in 
turn, with nonequilibrium situations. 

The authors appreciate the financial support of the 
SFB 458 and the University of Miinster as well as the 
partial support of DFG project MA 1508/8-1 and RFBR 
grants 06-01-04005, 05-01-00723, and 05-07-90248. 



2sin(^)r(gf^) 

v^ar(^) 



(14) 



follows directly form the comparison of the asymptotics 
of Px{Ax) determined by ([9]) and the asymptotics (fT2|) 
of the velocity extremum distribution. In particular, for 
a — 1.6 we have Ax ~ 1.6i?t in agreement with the 
simulation data. 

The developed model ([3|) actually gives us the imple- 
mentation of Levy flights at the "microscopic" level ad- 
mitting the notion of continuous trajectories, ndeed, fix- 
ing any small duration St of the Levy walker steps we can 
choose the time scale r of model ([3]) such that 6t ^ t 
and, as a result, receive the Levy statistics for the cor- 
responding spatial steps. Moreover, the found expres- 
sion ([9]) demonstrates the equivalence of all the systems 
in asymptotic behavior for which the parameters Va and 
r are related by the expression ?;"t"~^ = a. In some 
sense, all the details of the microscopic implementation 
of Levy flights are aggregated in two constants: the expo- 
nent a and the superdiffusion coefficient a. In particular, 
the characteristic scale of the walker displacement during 
time t is ~ (at)^/". 

Our approach has several immediate consequences. 
First of all, it yields an easily implementable procedure 
for the numerical simulation of Levy processes based on 
the simulation of the Langevin equations ([2]). Second, 
it seems to be possible to attack the yet unsolved prob- 
lem of the formulation of accurate boundary conditions 
for the generalized Fokker-Planck equations describing 
Levy processes in finite domains and heterogeneous me- 
dia. The crucial point of our treatment is the existence of 
quantities varying on three widely separated time scales 
(5t <C T ^ St. On time scales St the Langevin equation 
is updated. In the well-defined limit of small St the tra- 
jectory can be constructed with arbitrary precision. Fur- 
thermore, r is connected with the relaxation time of the 
variable v and sets the overall time scale of the model. Fi- 
nally, for St the variation of the position x is fully Marko- 
vian and the systems behaves according to the standard 
Levy flight scenario. A similar approach is the treatment 
of the Kramers-Fokker-Planck equation describing diffu- 
sion of particles, which is obtained from eq. (O for the 
case of purely additive noise g = const. The so-called 
Smoluchowski limit r — > leads to Einstein's diffusion 



[1] T.H. Solomon, E.R. Weeks, and H. L. Swinney, Phys. 

Rev. Lett. 71, 3975 (1993). 
[2] J.-P. Bouchaud, A. Georges, Phys. Rep. 195, 127 (1990). 
[3] G.M. Viswanathan, V. Afanasyev, S. Buldyrev, E. J. 

Murphy, P. A. Prince, H.E. Stanley, Nature 381, 413 

(1996). 

[4] D. Brockmann, L. Hufnagel, T. Geisel, Nature 439, 462 
(2006). 

[5] R.N. Mategna, H.E. Stanley, Nature 376, 46 (1995). 

[6] J. Klafter, M. F. Shlesinger, and G. Zumofen, Phys. To- 
day 49(2), 33 (1996). 

[7] Levy Flights and Related Topics in Physics, edited by 
M. Shlesinger, G. M. Zaslavsky, and U. Frisch (Springer 
Verlag, Berlin, 1994); Anomalous Diffusion, edited by A. 
Pekalski, K. Sznajd-Weron, Lecture notes in Physics Vol. 
519, (Springer Verlag. Berlin, 1999). 

[8] H. Risken, The Fokker-Planck Equation: Methods of So- 
lution and Applications 2-nd ed. (Springer- Verlag, Berlin, 
1996). 

[9] C. W. Gardiner, Handbook of Stochastic Methods for 

Physics, Chemistry, and the Natural Sciences, 3-rd ed. 

(Springer- Verlag, Berlin, 2004). 
[10] H. C. Fogedby, Phys. Rev. E 50, 1657 (1994). 
[11] A. Weron, K. Burnecki, S. Mercik, and K. Weron, Phys. 

Rev. E 71, 016113 (2005). 
[12] V. V. Yanovsky, A. V. Chechkin, D. Schertzer, and 

A. V. Tur, Physica A 282, 13 (2000). 
[13] D. Schertzer, M. Larcheveque, J. Duan, V. V. Yanovsky, 

and S. Lovejoy, J. Math. Phys. 42, 200 (2001). 
[14] I. M. Sokolov, Phys. Rev E 63, 011104 (2000). 
[15] H. C. Fogedby, Phys. Rev. Lett. 73, 2517 (1994). 
[16] J. Honkonen, Phys. Rev. E 53, 327 (1996). 
[17] D. Brockmann and T. Geisel, Phys. Rev. Lett. 90, 170601 

(2003). 

[18] D. Brockmann, L M. Sokolov, Chem. Phys. 284 409 
(2002). 

[19] S. Jespersen, R. Metzler, and H. C. Fogedby, Levy 
Flights in External Force Fields: Langevin and Frac- 
tional Fokker-Planck Equations, and their Solutions 
arXiv.org:cond-mat/9810176 (1998). 

[20] S. V. Buldyrev, S. Havlin, A. Ya. Kazakov, 
M. G. E. da Luz, E. P. Raposo, H. E. Stanley, 
and G. M. Viswanathanl, Phys. Rev. E 64, 041108 
(2001). 

[21] A. V. Chechkin, V. Yu. Gonchar, J. Klafter, and R. Met- 
zler, Europhys. Lett., 72, 348 (2005). 

[22] P. Imkeller and I. Pavlyukevich, J. Phys. A 39, L237 
(2006). 

[23] B. Dybiec, E. Gudowska-Nowak, and P. Hanggi, Phys. 



5 



Rev. E 73, 046104 (2006). 
[24] B. Dybiec, E. Gudowska-Nowak, P. Hanggi, Phys. Rev. 

E 75, 021109 (2007). 
[25] A. V. Chechkin, O. Yu. Sliusarenko, R. Metzler, J. 

Klafter, Phys. Rev. E 75, 041101 (2007). 
[26] S. Siegert, R. Friedrich, Phys. Rev. E 64, 041107 (2001). 
[27] G. Zumofen and J. Klafter, Phys. Rev. E 51, 2805 (1995). 
[28] P. M. Drysdale and P. A. Robinson, Phys. Rev. E 58, 

5382 (1998). 



[29] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 
[30] P. Levy, Theorie de I'Addition des Variables Aleatoires 

(Gauthier-Villars, Paris, 1954); B.V. Gnedenko and A.N. 

Kolmogorov, Limit Distributions for Sums of Random 

Variables, (Addison-Wesley, Reading, MA. 1954). 
[31] K. Lindenberg, K. E. Shuler, J. Freeman and T. J. Lie, 

J. Stat. Phys., 12, 217 (1975). 



Supplementary material to paper "Realizat ion of Levy flights as continuous processes" 



STOCHASTIC SYSTEM AND ITS GOVERNING 
EQUATION 

We consider continuous ID random walks governed 
by the following stochastic differential equations of the 
Hanggi-Klimontovich type [J, 0, 0| 



dx 
It 
dv 
'dt 



(g + l) 
2r 



(15) 
(16) 



where x is the position of a walker, v is its current ve- 
locity, ^(i) is white noise such that {^{t)^{t')) = 5{t — t'), 
and the function 



g{v) = \/vl+v^ 



(17) 



specifies the intensity of random Langevin forces. The di- 
mensionless coefficient a, the time scale r, and the char- 
acteristic velocity Va quantifying the additive component 
of the Langevin forces are the system parameters. The 
Levy flights arise when the coefficient a belongs to the 
interval 



1 < a < 2 



(18) 



which, thereby, is assumed to hold beforehand. 

For the given system the distribution function 
V {x — xo,v,VQ,t) obeys the following Fokker-Planck 
equation written in the kinetic form 



dV 
dt 



J_d_ 

2^dv 



dV 

g^{v)— + {a + l)vV 



^JW] (19) 



subject to the initial condition 
V {x — xo, V, vafi) — 6{x 



xo)5[v - vo) , (20) 



where, xo and are the initial position and velocity 
of the walker and, an addition, the system translation 
invariance with respect to the variable x is taken into 
account. 



VELOCITY DISTRIBUTION 
General relations 

It is the statistical properties of the walker velocity v 
that give rise to Levy flights. So the present section is 
devoted to them individually. The velocity distribution 
is given by the partial distribution function 



Py{v,VO,t) 



dxV {x — xo,v,vo,t) (21) 



which, by virtue of (|19p . obeys the reduced forward 
Fokker-Planck equation 



dt dv 



{v)—^ + {a + l)vP,, 
ov 



(22) 



written in the kinetic form whose right-hand side acts on 
the variable v. Simultaneously, the function Py{v,VQ,t) 
meets the backward Fokker-Planck equation 



dP, 2 d'p, dP, 



(23) 



written in the Ito form and acting on the variable vq (see, 
e-g-, 01 )■ The two equations are supplemented with the 
initial condition 



Py (w,wo,0) = 6{v - Vq) . 



(24) 



In particular, as stems from (j22p . the stationary velocity 
distribution -P^'(w) is the solution of the equation 



Qpst 

dv 



(a-Hl)wPf =0, 



(25) 



which together with the normalization condition 

dvP^\v) = l (26) 
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gives us the expression 



K 2 , 



0Fr(f) 



+1 



(27) 



where r(. . .) is the Gamma function. For the exponent a 
belonging to interval the first moment of the velocity 
V converges, whereas the second one diverges what actu- 
ally was the reason for specifying the region of a under 
consideration. 



First passage time problem and extremum 
distribution 



In order to establish some kinematic relationship be- 
tween the Levy type behavior exhibited by the given ran- 
dom walks on time scales t ^ t and properties of the 
velocity distribution we will make use of the first passage 
time statistics. The probability F{vo, d, t) for the walker 
with initial velocity wq such that \vo\ < "d to gain the ve- 
locity V = ±1^ for the first time at the moment t is directly 
described by the backward Fokker-Planck equation 
In particular its Laplace transform 



Jo 



obeys the equation (see, e.g., [1]) 



d^FL 



{a - l)vo 



dFr 



dvl ^" dvo 
subject to the boundary condition 



(28) 



(29) 



The introduced first passage time probability is necessary 
to analyze the extremum statistics. Namely we need the 
probability $(uo,^?,i) for the velocity pattern v{t) orig- 
inating from the point vq G +'&) to the get the ex- 
tremum equal to somewhen during the time interval 
t is related to the probability F{vo,-d,t) by the expres- 
sion 



d_ 



dt' F{vo,i},t') 



(30) 



or for the Laplace transforms 



1 d 

$l(wo,i?,s) = — —Fl{vo,^,s). 
s ov 



(31) 



To examine the characteristic properties of the first 
passage time statistics let us consider two limit cases, 
s ~^ and — > 00. Their analysis starts at the 
first step with the same procedure. Namely, we assume 
the function Fl(vq,'3, s) to be approximately constant, 
Fl{vo,i^, s) ~ Fo(i?, s) inside some neighborhood Qo of 
the origin vq — 0. For s — > it is the domain (—1?, i9) 
itself and Fo(i9, s) = 1 by virtue of ^Q^i. For -> 00 
the thickness of this neighborhood is much larger then 
Va as it will be seen below. Under such conditions equa- 
tion (|28p can be integrated directly inside the domain 
Qo with respect to the formal variable /(wo) — dF^/dyQ 
using the standard parameter-variation method. In this 
way taking into account that /(O) = due to the system 
symmetry we obtain the expression 



dvQ Va \vl J Jo (£.^ + 1) 

and for \vo\ ^ Va 

— V^rsFo i^,s) . (33) 



Expression ([33]) demonstrates us that, first, the imple- which can be rewritten as 
mentation of the limit case of small values of s (formally, 
s — > 0) is the validity of the inequality 

f„(^,,)»,„^^%M ^ rs^«l, ^l{s):=[^J Va»\vo\ (34) 

duo 
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or converting to the time dependence 



( - ) Va->\vo\ 



(35) 



So the characteristic velocity scale characterizing the first 
passage time probability and aggregating its time depen- 
dence is ^{t). As a consequence, the limit of small values 
of s is actually defined by the inequality \v\o <C ^{t). 
Correspondingly, the limit of large values of "i? is imple- 
mented by the inequality i) ^ d{t) or i? ^ ^l{s)- 

Second, for d ^ ^l{s) there is a region, namely, 
Va ^ \vo\ ^ ^{t) wherein the assumption FL{vo,'d, s) ~ 
Fq (d, s) holds whereas the derivative dF^ / dvQ scales with 
vq as dFj^/dvo oc |wo|" ^- This asymptotic behavior 
can be obtained also by analyzing the solution of equa- 
tion for I Vol » Va where g^{vo) — Vq. In this case 
equation admits two solutions of the form 



with 



~ a and 52 



91,2 



2ts 

a 



(36) 



The second solution is relevant to the function 
FlIvq, "i?, s) only within the crossover from Fl{vq, t?, s) oc 
Vq to Fl{vo,i}, s) ~ Fo{'d,s) and even in this region, i.e. 
ko| ^ ^{t) the derivative dF^/dvo is determined by its 
asymptotics Fi(uo,t9, s) oc Uq. For larger values of vq, 
i.e., |uo| ^ i?(i) the first passage time distribution is given 
by the expression 



Fl{vo,^, 



Nl 



(37) 



taking into account the boundary condition (P^)) . So we 
can write 



dFLivo,^,s) \vo\ 
~ a- 



dvo 



(38) 



also for |wo| < i?(t). Expressions ((33|) and p8|) de- 
scribe the same asymptotic behavior of the function 
Fl{vot'&,s). Thereby we can "glue" them together, ob- 
taining the expression for 



Fq{'&,s) 



' a+l 
, 2 



) 1 < 



0Fr (f ) Ts 



(39) 



It should be noted that this procedure is the kernel of 
the singular perturbation technique which will be also 
used below. Expression ((39)) immediately gives us the de- 
sired formula for the extremum distribution ^^(wo, 1?, s)- 
Namely, by virtue of (PT|) . for |wo| ^ "^Lis) we have 



$l(wo,i?,s) 



^r(^) 1 



TS- 



(40) 



and restoring the time dependence of the extremum dis- 
tribution from its Laplace transform the asymptotic be- 
havior for 3> ^{t) we get 



Hvo,^,t) 



'r(^)t v^ 



(41) 



Finalizing the present subsection we draw the conclu- 
sion that for I Wo I ^ 'i?(t) the extremum distribution is 
described by a certain function 



with the asymptotics 
$0 (0 = 



^r(H±i) 1 



(42) 



(43) 



Here the velocity scale i3(i) is given by expression (|35p . 
We remind that distribution (|4T|l describes the amplitude 
of the velocity extrema, so, as the velocity extrema them- 
selves are concerned their distribution is characterized by 
the function 



$(t;o,^,t) = ^$(«o,|^|,i) 



(44) 



because of the symmetry in the velocity fluctuations. 

It should be also noted that the asymptotics 
<I>(wo, ^9, t) oc for i3 > §(t) could be obtained im- 

mediately from equation (|28p . In fact, formally assuming 
I? — > 00 and taking into account the boundary condi- 
tion (|29p we can represent the solution of equation P5|) 
in form ([37)l for vq < d because, first, g^{vQ) = Wg in this 
case and, second, the function Fl^vo,!!), s) must be de- 
creasing with \vo\. It is the only one place where the vari- 
able d enters the function Fl{vo, t?, s), thus, for i9 ^ d{t) 



Fl{vo,i},s) oc 



1 



and relationship ([3T|l directly gives rise to 

1 



<i>{vo,^,t) oc 



GENERATING FUNCTION 

General relations 

To analyze the given stochastic process the generating 
function 

g{k,K,t) = <^exp I [rvk + {x - xq) k]^'^ (45) 

is introduced. As follows from the Fokker-Planck equa- 
tion pg]) it obeys the governing equation 
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subject to the initial condition 



g{k, >c, 0) = exp <^ — vok ^ . 



(47) 



At the origin fc = and k = Q function (US)) meets also 
the identity 



e(o,o,t)-i 



(48) 



which follows directly from the meaning of probability. 
In deriving equation (|46p the following relationships be- 
tween the operators acting in the spaces {x, v} and {>r, k} 

d i d i , d 

-E- ^ ^, TT ^ V ^ -^'"alTr 

ox VaT OV Vc OK 

as well as the commutation rule 

dk" "dk'^ 

have been used. 

The argument >f enters equation (|46p as a parameter; 
the given equation does not contain any differential oper- 
ator acting upon the function G{k, >c, t) via the argument 
K. This property enables us to pose a question about the 
spectrum of equation (j46p . where the variable x plays the 
role of a parameter. The desired eigenf unctions and their 
eigenvalues 



{Vl/A(fc|x)}, {k{>,)} 

obey the equation 



(49) 



-2A*A — f A: 



dk V dk 



dk 



(50) 

In deriving equation (jSOp the time dependence 
exp(— Af/r) corresponding to eigenfunctions has 
been assumed. 

In these terms the solution of equation (|46|) is reduced 
to the series 

^;(fc, X, i) = ^ /a {x, |f o) *a {k\x) exp I - A (k) ^ | 

(51) 

whose the coefficients {f {x,\vo)} meet the equality 

^/A(>f, Ivo)^'A(fck) =expi — wofcl (52) 



steaming from the initial condition (j47[) . In agreement 
with the results to be obtained, the spectrum of the 
Fokker-Planck equation (|46|) is bounded from below by a 
nondegenerate minimal eigenvalue Amin (><) > whereas 
the other eigenvalues are separated from it by a final gap 
of order unity. 



g{k, K, t) = /mill {x, \vo) ^fmin {k\x) CXp \ -A nun (^) 



So, as time goes on and the inequality t t holds, 
the term corresponding to the minimal eigenvalue will 
be dominant and sum (1511) is reduced to 



t 



(53) 

on large time scales. Here ^min {k\>c) is the eigenfunction 
of the eigenvalue Amin (>?) • 

Whence several consequences follow. First, the iden- 
tity holds at any time moment, thereby 



Amin (0) = . 



(54) 



Second, in the limit case t ^ t the system has to "forget" 
the value vq of initial velocity, so the coefficient /min (>^) 
does not depend on the argument vq and, therefore, can 
be aggregated into the function ^'min In this way 

the initial condition expansion (I52|) reads 



*mm(fck)+ 



X — ^ I ^ 

/a (x, |t^o) *A (fc|>f) = exp <^ — vak 



A>A „ 



(55) 
for any vq. 

The terms in sum (|55p with A > Amin determine 
the dependence of the generating function Q{k^>c,t) on 
the initial velocity wq, so, the corresponding coefficients 
/a {x, \ v{)) must depend on wq. Finding the first deriva- 
tive of both the sides of this equality with respect to vq 
we have 



tt-Za (>^, bo) *A (fc|>^) = — /cexpj — Vok 



Whence it follows that, third, the eigenfunctions 
^Pa {k\x) for A > Amin must exhibit the asymptotic be- 
havior vPa (fcjx:) ^ as A: ^ because of their linear 
independence. Fourth, setting fc = in expression ([55|) 
we get the conclusion that the eigenfunction 5'min {k\x) 
has to meet the normalization condition 



{Q\x)^\ at A: = 0. 



(56) 



Summarizing the aforementioned we see that on large 
time scales t ^ t the desired asymptotic behavior of the 
given system is described by the generating function 

g{k, X, t) = *min {k\x) exp I - Amin (^) ^ | ' 



(58) 



and by virtue of ([56 



g{0, X, t) = exp \ -Amin {x) - 



In what follows the calculation of the eigenvalue Amin (x) 
will be the main goal. 
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The given random walks should exhibit the Levy flight 
behavior on large spatial and temporal scales, i.e. x 3> 
VaT and T. It allows us to confine our analysis to the 
limit of small values of x, i.e. assume that \>c\ <C 1, where 
also the eigenvalue Amin(><) ^ 1. In this case the spec- 
trum of equation may be studied using perturbation 
technique, where the term 



plays the role of perturbation. 



dk 



(59) 



nonzero value at the origin fc = and, thereby, enables 
the eigenfunction problem (|63p to be considered within 
i^-space. 

Second, for /32 = a/2 equation ([50)1 is reduced to the 
modified Bessel differential equation 



9a. 2 



dk 



2A 



<^A.2-0. (65) 



Since the desired eigenfunctions should decrease as A: — *■ 
00 the solution of equation ((GS]) is given by the modified 
Bessel function of the second kind 



Zero-th approximation. Spectral properties of the 
velocity distribution 



The zero-th approximation of (j46p in perturbation (I59p 
matches the case >f = 0, where the generating func- 
tion (|45)l actually describes the velocity distribution ((2T|) . 
Setting X — reduces the eigenvalue equation ([50]) to the 
following 



2A$> 



d 
dk 



dk 



{a + l)k 



dk 



where 



*A(fc) = *A (A:|0) and A = A (0) 



(60) 



(61) 



Having in mind different goals we consider the conversion 
of equation ()60p under the replacement 



$A(fc) - \kf' </)A,.(fc) 



for two values of the exponent (3i. 

First, for /?i = {a + l)/2 equation ([60]) is converted 
into 



2A(/)A,i = 



d 

dk 



k' 



dk 



1 



(63) 

The operator on the right-hand side of equation ([55]) is 
Hermitian within the standard definition of scalar prod- 
uct. So all the eigenvalues A are real numbers and the 
corresponding eigenfunctions form a basis. It should 
be noted that the given conclusion coincides with the 
well known property of the Fokker-Planck equations with 
the detailed balance @. In addition the eigenfunctions 
0A,i (k) can be chosen so that the identity 



dk(f>l^^ (k) (f>y,i (k) = Sxx' 



(64) 



holds for all of them except for the eigenfunction 
4>min (k) corresponding to the minimal eigenvalue Amin = 
Amin (0) = by virtue of ([54]). We note that the latter 
eigenfunction describes the stationary velocity distribu- 
tion ()27p and its normalization is determined by condi- 
tion (|56p . Treating the eigenfunction ^min (k) individu- 
ally releases the remainders from the necessity to take a 



^xAk) 'X K,{\k\) 



(66) 



with the order v 



KA\k\) 



4" 



2A because 



2\k\ 



as 



oo 



for any value of the parameter u 

Whence it follows that there are no eigenfunctions with 
A < ^a^ and A 7^ 0. Indeed, when A < the function 



$(fc) |fcp"if,(|fc|)cx |fc| 



for \k\ < 1 



diverges as fc — > 0. In the region < A < ^a'^ the 
corresponding eigenfunctions 



Mk) = ^xAk)\k 



cx \k\-^K,i\k\) 



(62) 

would give rise to a strong divergency in the normaliza- 
tion condition (|64p. When A > the solution of equa- 



tion ((65)) is described by the modified Bessel functions of 
pure imaginary order which exhibit strongly oscillatory 
behavior as x: ^ and describe the continuous spec- 
trum of the Fokker-Planck equation (j46|l for >«: = 0. Due 
to result the eigenfunction <i>i„in {k) corresponding 
to the eigenvalue A = and meeting the normalization 
condition (l56t is of the form 



*min (fc) 



^(^fc^i^tdfcl) 



' 2 j r(^) 



O(fc^). (67) 



In deriving expression (j67p the following expansion of the 
function K^{k) has been used 



K.{\k\)^ 



ifcr 



r(i + z.) 



Oik' 



(68) 

which is justified for the order < < 1 (see, e.g., 
Ref. 01 )■ The latter inequality holds due to the adopted 
assumption (|18p about the possible values of the param- 
eter a. 
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Expression ((67)) finalizes the analysis of the zero-th ap- 
proximation. Summarizing the aforementioned we draw 
the conclusion that at x = the spectrum of the Fokker- 
Planck equation ([^5]) for the generating function ([IS]) 
does contain zero eigenvalue Amin (0) = correspond- 
ing to eigenfunction (j67p which is separated from higher 
eigenvalues by a gap equal to (in units of t). Wc 

note that the given statement is in agreement with the 
conclusion about the spectrum properties for a similar 
stochastic process with multiplicative noise [1, [13] ■ 



side of equation ((50)) is also ignorable in comparison with 
its left-hand side. Under these conditions the eigenvalue 
equation (|50)) is reduced to the following 



2A, 



dk 



Ik 



(69) 



for the function ip {k\>c) = 1 — ^min (fc|>f)- In the given 
case the singular perturbation technique is implemented 
within the replacement k = Qx converting equation 
into one of the form 



The eigenvalue Amin(x) for |x| ^ 1. Singular 
perturbation technique 



2A„ 



When >ir 7^ the perturbation term ((59)) mixes the 
eigenfunctions of zero-th approximation and, as a re- 
sult, the eigenfunctions ^x{k) with A > contribute 
also to the eigenfunction ^'min However, because 

their eigenvalues are about unity or larger, A > 1, 
the perturbation can be significant only in the domain 
1^1 ^ 1^1- Outside this domain the perturbation is not 
essential and the eigenfunction ^'min (^|^) practically co- 
incides with its its zero-th approximation $min(^)- So 
in the case when |>f| ^ 1 there should be an inter- 
val \k\ <^ \k\ <^ 1 where, on one hand, the eigenfunc- 
tion 5'niin {k\x) can be already approximated by <i>inin (k) 
and, on the other hand, the expansion ((67)) still holds, 
in particular, ^'min(fc|>£') ft; 1 in this region. Leaping 



ahead, we note that Aj, 



I" so inside the subinter- 



val \x\ <^ \k\ <^ \h\^ the last term on the right-hand 



subject to the effective "boundary" conditions by virtue 
of ((63 



C^±oo. (71) 



In some sense the condition ((7T)) "glues" the asymptotic 
behavior of the eigenfunction ^!^i^{k\x) resulting from 
its properties for sufficiently large values of k together 
with one stemming from small values of fc, in this case, 
specified by the solution of equation ((TT)) . Exactly such 
a procedure is the essence of the singular perturbation 
technique. 

Equation ()70p with respect to the variable dip/d( can 
be solved directly using the standard parameter- variation 
method. In this way we get for C < 



d4 ^ ,.,a-l 

dc 



exp [ - I C-c 



^ A 



r^^expi ^-ndc 



2/C 



(72) 



and for C > 



dip i^ia-l / 2 

dc = U 



a 



+ OC 



2/C 

2^-"A^inJ r"'exp(-e)de 



(73) 



where the constants C±oo specify the asymptotic behav- 
ior of the derivative 

dip iz-ia-i^ ^ , 

~ ICI C:±oo as C ^ ±oo 

and according to condition ()7ip 



Expression ((75)) diverges as C ^ unless the equality 



- A,„in / exp(-Or"'rfe = 



r(i-^) 



holds, whence we find the desired expression for the 
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eigenvalue A„ 



-^min 



r(a)r(f) 



(75) 



Expression ([75)l finalizes the analysis of the generat- 
ing function (|45p . In particular, together with expres- 
sion ([58|l it gives the desired formula for the generating 
function 



exp 



{x — Xq) k 



VaT 



exp • 



r(a)r(f) 



(76) 



spect to the walker displacement x — Xq and its velocity 
extrema we get the relationship between their charac- 
teristic scales 



X{t) = i){t)T 



(79) 



and the asymptotic equivalence within the replacement 
{x — Xq) = dT, where 



T ■ 



2sin(^)r 



\ 2 ; 



0Far(^) 



(80) 



demonstrating the fact that the given random walks ex- 
hibit Levy flight statistics on time scales t ^ t. Ex- 
pression (I76p in turn gives us the asymptotics of the x- 
distribution function 

Px{x - xo,vo,t) ^ / dvV {x - xo,v,vo,t) 



for |a; — a;o| 3> x(t) in the form 
sm{^) ar( 



Pxix - xo,t) = 
where the length 



2~a ^ 
2 J 



c"(i) 



x{t) = {at)' 



\x - xo\ 



W (77) 



(78) 



with cr = w"t" ^ specifies the characteristic scales of the 
walker displacement during the time interval t. 



THE LEVY FLIGHT BEHAVIOR AND THE 
EXTREMUM STATISTICS OF THE WALKER 
VELOCITIES 

Comparing expressions ([H]) and ([77]) describing the 
asymptotic behavior of the given random walks with re- 



The obtained expressions allow us to consider the long 
distance displacements of the walker within the time in- 
terval t to be implemented during one spike of duration 
r in the pattern v{t) that has the maximal amplitude. In 
particular, for a = 1.6 the ratio T/r ~ 1.6 
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